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Abstract. We discuss quantum graphs consisting of a compact part and semiinfinitc 
leads. Such a system may have embedded eigenvalues if some edge lengths in the 
compact part are rationally related. If such a relation is perturbed these eigenvalues 
may turn into resonances; we analyze this effect both generally and in simple examples. 



Perturbed quantum graphs with rationally related edges 



2 



1. Introduction 

Quantum graphs have attracted a lot of attention recently. The reason is not only 
that they represent a suitable model for various microstructures, being thus of a direct 
practical value, but also that they are an excellent laboratory to study a variety of 
quantum effects. This comes from a combination of two features. On one hand these 
models are mathematically accessible since the objects involved are ordinary differential 
operators. On the other hand graphs may exhibit a rich geometrical and topological 
structure which influences behaviour of quantum particle for which such a graph is 
a configuration space. There is nowadays a huge literature on quantum graphs and, 
instead of presenting a long list of references, we restrict ourselves to mentioning the 
review papers |Ku04,05 lKu08j as a guide to further reading. 



One important property of quantum graphs is that - - in contrast to usual 
Schrodinger operators — the unique continuation principle is in general not valid for 
them: they can exhibit eigenvalues with compactly supported eigenfunctions even if 
the graph extends to infinity. This property is closely connected with the fact that 
eigenvalues embedded in the continuous spectrum are on quantum graphs by far less 
exceptional then for usual Schrodinger operators. A typical situation when this happens 
is when the graph contains a loop consisting of edges with rationally related lengths 
and the eigenfuction has zeros at the corresponding vertices, which prevents it from 
"communicating" with the rest of the graph. 

On the other hand, since such an effect leans on rational relations between the edge 
lengths, it is unstable with respect to perturbations which change these ratios. The 
resolvent poles associated with the embedded eigenvalues do not disappear under such 
a geometric perturbation, though, and one can naturally expect that they move into 
the second sheet of the complex energy surface producing resonances. The aim of the 
present paper is to discuss this effect in a reasonably general setting. 

We consider a graph consisting of a compact "inner" part to which a finite number of 
semiifinite leads are attached. We assume a completely general coupling of wavefunctions 
at the graph vertices consistent with the self-adjointness requirement. As a preliminary 
we will show, generalizing the result of |EL06j, that we can speak about resonances 
without further adjectives because the resolvent and scattering resonances coincide in 
the present case. We also show how the problem can be rephrased on the compact graph 
part only by introducing an effective, energy-dependent coupling. 

After that we formulate general conditions under which such a quantum graph 
possesses embedded eigenvalues in terms of the graph geometry (edge lengths) and the 
matrix of coupling parameters. The discussion of the behaviour of embedded eigenvalues 
is opened by a detailed analysis of two simple examples, those of a "loop" and a "cross" 
resonator graphs. Here we can analyze not only the effect of small length perturbations 
but also, using numerical solutions, to find the global pole behaviour and to illustrate 
several different types of it. Returning to the general analysis in the closing section, 
we will derive conditions under which the eigenvalues remain embedded, and show that 
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"nothing is lost at the perturbation" in the sense that the number of poles, multiplicity 
taken into account, is preserved. 

2. Preliminaries 

2.1. A universal setting for graphs with leads 

Let us consider a graph V consisting of a set of vertices V = {Xj : j G /}, a set of finite 
edges £ = {£j n '■ (Xj, X n ) G Jc C fx /} and a set of infinite edges = {£joo : Xj G Ic} 
attached to them. We regard it as a configuration space of a quantum system with the 
Hilbert space 

^ = £ 2 ([0,y)© £ 2 ([0,oo)). 

the elements of which can be written as columns ip = (fj : Lj G C, gj : £joo £ £oo) T - 
We consider the dynamics governed by a Hamiltonian which acts as — d 2 /dx 2 on each 
link. In order to make it a self-adjoint operator, boundary conditions 

{Us-Wj + nUj + 1)^ = (i) 

with unitary matrices Uj have to be imposed at the vertices Xj, where tyj and ^ are 
vectors of the functional values and of the (outward) derivatives at the particular vertex, 
respectively. In other words, the domain of the Hamiltonian consists of all functions 
in W 2 ' 2 (C © Coo) which satisfy the conditions ([TJ. We will speak about the described 
structure as of a quantum graph and as long as there is no danger of misunderstanding 
we will use for simplicity the symbol T again. 

While the model is simple dealing with a complicated graph may be nevertheless 
cumbersome. To make it easier we will employ a trick mentioned to our knowledge for 
the first time in |Ku08] passing to a graph To in which all edge ends meet in a single 
vertex as sketched in Fig. [T] the actual topology of T will be then encoded into the 
matrix which describes the coupling in the vertex. 

To be more specific, suppose that T described above has an adjacency matrix CV,- 
and that matrices Uj describe the coupling between vectors of functional values ^ and 
derivatives ^ at Xj. This will correspond to the "flower- like" graph with one vertex, 
the set of loops isomorphic to C and the set of semiinfinite links which does not 
change; coupling at the only vertex of this graph is given by a "big" unitary matrix U. 

Denoting iV = card L and M = card C^, we introduce the (2iV + M)-dimensional 
vector of functional values by $ = , . . . , \l/^ ai . dV ) T and similarly the vector of 
derivatives ^' at the vertex. The valency of this vertex is M + J2i.j Cij — 2iV + M. 
One can easily check that the conditions (fij) can be rewritten on T using one 
(2N + M) x (2iV + M) unitary block diagonal matrix U consisting of blocks Uj as 

{U - 1)^! + i{U + 1)^' = 0; (2) 

the equation (l2l) obviously decouples into the set of equations (fl| for and \E^. 
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Figure 1. The model To for a quantum graph T with N internal finite edges and M 
external links 



Since neither the edge lengths and the corresponding Hilbert spaces nor the operator 
action on them are affected and the only change is a possible edge renumbering the 
quantum graph Tq is related to the original T by the natural unitary equivalence and 
the spectral properties we are interested in are not affected by the model modification. 



2.2. Equivalence of the scattering and resolvent resonances 

As another preliminary we need a few facts about resonances on quantum graphs. In 
[EL06] we studied the situation where to each vertex of a compact graph at most one 
external semi-infinite link is attached; we have demonstrated that the resonances may 
be equivalently understood as poles of the analytically continued resolvent, (H — A id) -1 , 
or of the on-shell scattering matrix. Here we extend the result to all quantum graphs 
with finite number of edges, both finite and semi-infinite: we will show that the resolvent 
and scattering resonances again coincide. The above described "flower-like" graph model 
allows us to give an elegant proof of this claim. 

Let us begin with the resolvent resonances. As in [E L06j the idea is to employ an 
exterior complex scaling; this seminal idea can be traced back to the work J.-M. Combes 
and coauthors, cf. |AC71j . and its use in the graph setting is particularly simple. Looking 
for complex eigenvalues of the scaled operator we do not change the compact-graph part: 
using the Ansatz fj (x) = dj sin kx + bj cos kx on the internal edges we obtain 

/j(0) = bj, fj(lj) = a j sm klj + bj cos klj, (3) 

fj(Q) = kctj, — f'jilj) = —kaj cos klj + kbj sin klj, (4) 
hence we have 

(5) 
(6) 

On the other hand, the functions on the semi- infinite edges are scaled by gje(x) = 
e e / 2 gj(xe 9 ) with an imaginary 9 rotating the essential spectrum of the transformed 
(non-selfadjoint) Hamiltonian into the lower complex halfplane so that the poles of the 
resolvent on the second sheet become "uncovered" if the rotation angle is large enough. 
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The argument is standard, both generally and in the graph setting |EL06] . so we skip 
the details. In particular, the "exterior" boundary values are given by 

g j (0) = *-° / *9i*> <£(0) = ifce-^W (7) 

Now we substitute eqs. @, @ and into @. We rearrange the terms in \1/ 
and in such a way that the functional values corresponding to the two ends of each 
edge are neighbouring, and the entries of the matrix U are rearranged accordingly. This 
yields 



{U-I)C x {k) 



a i 
bi 

«2 



on 
e^gw 



V e-^g Me J 



+ ik{U + I)C 2 (k) 



a i 
bi 

«2 



on 
e- e ' 2 gie 



\ e- 9/2 9M e J 



(8) 



f\k),c[ 2) ( 
respectively, where 



where the matrices Ci, C 2 are given by Ci(fc) = diag (C[ L) (k), CY\k), . . . , (A;), Imxm) 
and C 2 = diag (C^\k),C^\k), . . . , C^(A:), z/mxm), 

tW/^ = ( 1 ^ 

sin feL cos fcL / 



C[ 3 \k) 



Ci 3 \k) 



cos fcL- 





sin kh 



and Imxm is a M x M unit matrix. 

The solvability condition of the system ^ determines the eigenvalues of scaled 
non-self adjoint operator, and mutatis mutandis the poles of the analytically continued 
resolvent of the original graph Hamiltonian. 

The other standard approach to resonances is to study poles of the on-shell 
scattering matrix, again in the lower complex halfplane. In our particular case we 
choose a combination of two planar waves, gj = cje~ tkx + dje tkx , as an Ansatz on the 
external edges; we ask about poles of the matrix S = S(k) which maps the vector of 
amplitudes of the incoming waves c = {c n } into the vector of the amplitudes of the 
outgoing waves d = {d n } by d = Sc. The condition for the scattering resonances is then 
det S^ 1 = for appropriate complex values of k. The functional values and derivatives 
at the vertices are now given by 



= c 5 + d v 9j(Q) = i k (dj - cj). 
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together with eqs. ([3])-(|4]). After substituting into ^ one arrives at the condition 



{U-I)C^k) 



\ 



a i 
h 

0,2 

b N 
ci + d\ 



\ c M + d M ) 



+ ik{U + I)C 2 (k) 



a. i 

h 

b N 
d\ — c\ 



\ d M - c M / 



0. 



Since we are interested in zeros of det S \ we regard the previous relation as an equation 
for variables a,j, bj and dj while Cj are just parameters, in other words 



[{U-I)C x {k) + ik{U + I)C 2 {k)\ 



( CLi \ 

a 2 

b N 
d\ 

\d M ) 



-{U-^C^+ikiU + I^ik)} 



( \ 




Cl 

V cm j 



Eliminating the variables aj, bj one can derive from here a system of M equations 
expressing the map S ,_1 c/ = c. The condition under which the previous system is not 
solvable, what is equal to det 5* _1 = 0, reads 

det \{U - I) d(k) + ik{U + 7) C 2 (k)} = (9) 

being the same as the condition of solvability of the system ([8]); this means that the 
families of resonances determined in the two ways coincide. 



2. 3. Effective coupling on the finite graph 

The study of resonances can be further simplified by reducing it to a problem on the 
compact subgraph only. The idea is to replace the coupling at the vertex where external 
semi-infinite edges are attached by an effective one obtained by eliminating the external 
variables. Substituting from ^ into eqs. ^ we get 

f h \ ( A \ 



(U-I) 



f2N 

e -fl/2 



9w 



\ e~ e/2 9Me ) 



+ (U + I) diag(i, ... ,i, -k, . . . , -k) 



/2JV 
e-0/2 



9w 



0. (10) 
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We consider now U as a matrix consisting of four blocks, U = ( 1 2 ] , where U\ 

is the 2N x 2N square matrix referring to the compact subgraph, £7 4 is the M x M 
square matrix related to the exterior part, and U 2 and U3 are rectangular matrices of 
the size M x 2N and 2N x M, respectively, connecting the two. Then the previous set 
of equations can be written as 

V(A, ...,f 2N ,f f 2N) e' e ' 2 g ie , . . . , e- e l 2 g M6 f = 0, 

where 



V 



Ux-I z(C/i + /) (1 - k)U 2 \ 

v tf 3 iU 3 (l-k)U 4 -(k + l)I ) 

If the matrix [(1 — fc)?7 4 — (k + 1)] is regular, one obtains from here 



(e ' gxe,..., e ' 



)T = -[(!_ A;)f/ 4 - (k + m^Uzih + i/i, ...,f 2N + if 2N ) r 



and substituting it further into (10) we find that the following expression, 

[u x - 1 - (1 - k)u 2 [(i - k)u 4 - (k + 1)/]" 1 ^} (A, • • • , / 2 iv) T + 

+z + J — (1 — fc)t/ 2 [(l - A;)£/ 4 - (fc + l)/]- 1 ^} (/!,-•■, /^) T = • 

must vanish. Consequently, elimination of the external part leads to an effective coupling 
on the compact part of the graph expressed by the condition 

(U(k) - /)(/!, . . . , f 2N ) T + i(U(k) + /)(/{, . . . , f 2N ) T = , 

where the corresponding coupling matrix 

U(k) = C7i - (1 - k)U 2 [{\ - k)U 4 — (k + l)/]" 1 ^ (11) 
is obviously energy-dependent and, in general, may not be unitary. 

3. Embedded eigenvalues for graphs with rationally related edges 

As mentioned in the introduction, quantum graphs of the type we consider here have 
the positive halfline as the essential spectrum, and they may have eigenvalues with 
compactly supported eigenfunctions embedded in it. 

3.1. A general result 

We will focus on graphs which contain several internal edges of lengths equal to integer 
multiples of a fixed Iq > 0. In the spirit of the previous section we restrict ourselves only 
to compact graphs remembering that the presence of an exterior part can be rephrased 
through an effective energy-dependent coupling replacing the original U by the matrix 
U(k) defined above. 



Following Sec. 2.1 we model a given compact T by T having only one vertex and 



iV finite edges emanating from this vertex and ending at it. The coupling between the 
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edges is described by a 2N x 2N unitary matrix U and condition Suppose that 
the lengths of the first n edges are integer multiples of a positive real number l Q . Our 
aim is to find out for which matrices U the spectrum of the corresponding Hamiltonian 
H = Hu contains the eigenvalues k = 2mir/lo, m6N. 

Since our graph is not directed it is convenient to work in a setting invariant with 
respect to interchange of the edge ends. To this aim we choose the Ansatz 

tyj (x) = Aj sin k(x- lj 1 2) + Bj cos k(x - lj /2) . 



on the j-th edge. Subsequently, one gets 




sin 



cos 



cos 



cos 



cos 



sin 



sin 




The eigenvalue condition, expressed in terms of solvability of the system (|2]), is given by 
det[UD 1 (k) + D 2 (k)] = 0, (12) 



where 



Di(k) = 



/kl-i . -i kl-i 
— sin + ik cos 

sin — %k cos 



\ 



o 





D 2 (k) 



( 



\ 



sin ^ + ik cos 
- sin ^1 - ik cos sh. 



cos + ik sin 

fcTi i -i - fcTi 
cos + ik sm 






- cos + sm 

fcti i -7 ■ kl i 

- cos + ife sm ^ k 



fcZ AJ , ■ 7 kl AT 

sin + ik cos — 

kl ]\r ■, kl r\r 

sin —g*- — %k cos 



\ 



kl i\r . ■ j kl nj 

sin — + z/c cos — y~ 

kl i\r -j kl ' 

— «fc 



fcZ M i W - fcZ ?V 

cos + r/e sin — 
cos % + *fc sin / 



— COS 



— sin 



COS 



klw 
- COS — ^- 



■ i ■ kl /v 

-rj r ««: sin -j ?L 

-f - + ifcsin / 



For a future purpose, let us rewrite the spectral condition (12) in the form 



det(C(fc) + S(k)) = 0, where the matrix C(k) contains terms with cos -y- and S(k) 
contains those with sin k ^-. Hence all the entries in the first 2n columns of S(k) vanish 
for k = 2m7i/lo, m G N while the others can be nontrivial. Similarly, all the entries in 
the first 2n columns of C(k) are for k = (2m + 1)tt/1 , m G N equal to zero. The entries 
of the "cosine" matrix are 

kl ■ kU 

Ci,2j-i(k) = (u ii2 j-i - u it 2j)ik cos + (5 ij2 j-i - 5 it2 j) " 



S^2j)ik cos 

Ci,2j{k) = {Ui t2 j-l + U i)2 j) cos ^ - (S i)2 j-1 + S ij2j ) COS 



klri 



kl q 



Similarly, the entries of S(k) are 

Si,2j-l(k) = 



klj 
2 



u it2 j-i + Ui,2j) sin — - + {S ij2 j~i - S ij2j ) sin 



kl q 



2 ' 

kl- kl- 

Si,2j(k) = (iti,2i-i + u it2j )iksm — 1 + (f^-i + 5^ 2j )iksm . 



First of all, let us consider the situation when smklo/2 = 0. 
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Theorem 3.1. Let a graph Tq consist of a single vertex and N finite edges emanating 
from this vertex and ending at it, and suppose that the coupling between the edges is 
described by a 2N x 2N unitary matrix U and condition Let the lengths of the first 
n edges be integer multiples of a positive real number Iq. If the rectangular 2N x 2n 
matrix 



( 



M, 





U12 ~ 1 


U13 


u 14 




-1 


Ul )2n 


U 2 l ~ 1 


U22 


U23 


U24 


U2,2n~ 


-1 


U 2 ,2n 


U31 


U32 


U33 


U34 - 1 ■ ■ 


^3,2n- 


-1 


U 3 ,2n 


«41 


«42 


M43 - 1 


W44 


^4,2n- 


-1 


u 4,2n 



V 



U2N-l,2n-l 
U2N.2n-l 



U 2 N-l,2n 
U2N,2n 



(13) 



U 2 N-1,1 U 2 N-1,2 U 2 N~l,3 U2N-1,4 
U2N,1 U2N,2 U2N,3 U-2N,4 

has rank smaller than 2n then the spectrum of the corresponding Hamiltonian H = Hjj 
contains eigenvalues of the form e = 4m 2 ir 2 /ll with m e N and the multiplicity of these 
eigenvalues is at least the difference between 2n and the rank of M even . 



Proof. The condition (12) is clearly satisfied if the rectangular matrix containing only 
the first 2n columns has rank smaller than 2n, because then some of the columns of 
matrix C(k) + S(k) are linearly dependent. Since all the entries of the first 2n columns 
of S(k) contain the term sin klj/2, which disappear for kl = 2m7i, one can consider the 
matrix C(k) only. Dividing some of the columns of C(k) by appropriate nonzero terms, 
which is possible since cosklj/2 7^ for sinfc/ /2 = 0, and subtracting them from each 
other does not change the rank of the matrix. This argument shows that the rank of 
matrix M even must be smaller than 2n in order to yield a solution of the condition (12) 
and that the multiplicity is given by the difference. □ 

It is important to notice that the unitarity of U played no role in the argument, 
and consequently, one can obtain in this way embedded eigenvalues e = 4m 2 7i 2 /ll for a 
graph containing external links, however, the matrix M even (k) defined in analogy with 



(13) must have rank smaller than 2n for all values of k. 

Mathematically speaking the described case does not involve only cases where the 
original graph V contains a loop with rational rate of the lengths of the edges. Choosing 
appropriate U one can find such eigenvalues also for graphs where the edges of F with 
lengths equal to integer multiples of l are not adjacent. This corresponds, however, 
to couplings allowing the particle to "hop" between different vertices which is not so 
interesting from the point of view of the underlying physical model. 

A similar claim can be made also for Hq equal to odd multiples of 7r. 
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Figure 2. A loop of the edges with rational rate of their lengths 



Theorem 3.2. If under the same assumptions as above, the rectangular 2N x 2n matrix 



( 



odd 



Mil 


U12 + 1 


U13 


Uu 


Ul,2n- 


-1 


«l,2n 


U 2 1 + 1 


U22 


U23 


U24 


U2,2n- 


-1 


u 2,2n 


«31 


U 3 2 


U33 


U34 + 1 ■ • 


U 3 ,2n- 


-1 


U3,2n 


W 4 1 


UA2 


M43 + 1 


U44 


U4,2n- 


-1 


U4,2n 



\ 



(14) 



U2N-1,1 U 2 N-1,2 U 2 N~1,3 U 2 N-1,4: ' ' ' U 2 N-l,2n-l U2N-l„2n 
\ U2NA U2N,2 U-2N,3 U2N,4 ' ' ' U2N,2n-l U2N,2n J 

has rank smaller than 2n then the spectrum of the corresponding Hamiltonian H = Hjj 
contains eigenvalues of the form e = (2m + 1) 2 7t 2 //q with m 6 N and the multiplicity of 
these eigenvalues is at least the difference between 2n and rank of M a &&. 

We skip the proof which is similar to the previous one, the change being that the 
roles of the matrices S(k) and C(k) are interchanged. We also notice that similarly as 
above the results extends to graphs with semi-infinite external edges. 



3.2. A loop with 5 or 5' s couplings 

As mentioned above a prime example of embedded eigenvalues in the considered class of 
quantum graphs concerns the situation when T contains a subgraph in the form of a loop 
of n edges with the lengths equal to integer multiples of Iq. We denote by Uj, j — 1, . . . n, 
the unitary matrices describing the coupling at the vertices of such a loop and by U n+ i 
the unitary matrix which describes the coupling at all the other vertices of the graph 
- cf. Fig. [2j The unitary matrix which describes the coupling on the whole graph, in 
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the sense explained in Sec. 2.1, is 



U 









u 2 



V o o 



\ 


u n+1 J 



Let us further restrict our attention to the case when the coupling in the loop vertices 
is invariant with respect to the permutation of edges, i.e. suppose that matrices Ui, 
. . . , U n can be written as Uj = a,jj + bjl, where J is a unit matrix, J is a matrix 
with all entries equal to one and aj and bj are complex numbers satisfying \bj\ = 1 and 
\bj + ajdeg Xj\ = 1 to make the operator self-adjoint — cf. |ET07j . 



Recall that in order to use Theorems 3.1 and 3.2 one has to rearrange the columns 
and rows of the unitary matrix U accordingly. The first 2n entries in \l/ and 
correspond to the edges with rational rates of their lengths. Therefore, appropriate 
permutations of columns and rows of U must be performed: the first two columns 
should correspond to the first edge of the loop (from the vertex 1 to 2), the second two 

' U T \ 



columns to the second edge, etc. The rearranged coupling matrix is thus 
with 



a 



n+l 



+ &i 

a 2 + i>2 a.2 
a 2 a 2 + b 2 





oi 
ai 









02 





CI 2 

a 2 



0,1 ax 





+ b n a 




+ £>„ 
ai + 6i ai 
oi ai + bi 



0,1 









ai 
a 2 
a 2 







ai 

ai 



ai+bi 
02 + b 2 



a 2 





"2 
"2 





02 



02 + 6 2 



V 

The corresponding matrix M even is constructed in the way described in the previous 
section. It consists of a nontrivial 2n x 2n part (U T with added -l's) and degXi — 2 
copies of the row (a 1; 0, . . . , 0, ai), deg^ 2 — 2 copies of the row (0, a 2 , a 2 , . . . , 0), etc., 
and, finally, its last degX n+ i rows have all the entries equal to zero, hence the total 
number of its rows is 2N as required. 
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If all the Oj's are nonzero, the condition rankM e 



< In simplifies to 



rank 



bl 


-1 





■ 


■ 








-1 







■ 


• 














b 2 


-1 • 


• 








n 
u 


n 
u 


— 1 


h 
3 


n 
u 


n 
u 


n 
u 











■ 


• b n 

















■ 


■ 


K 


-1 











■ 


■ 


-1 




1 








■ 


■ 





1 





1 


1 


■ 


■ 

















■ 


■ 1 


1 






< 2n. 



It is easy to see that the assumptions of Theorem 3.1 giving rise to eigenvalues 
corresponding klo = 27cm are satisfied in the case bj = —1, Vj G {1,...N}, which 
corresponds to 5-couplings. The counterpart case, bj = 1, corresponding to 5' s couplings 
leads to the requirement 



rank 



/ 1 1 

1-10 

110 

1-1 



V 





-10 



\ 









1 -1 
1 1 
1 



< 2n 



which is satisfied if and only if the number of the edges in the loop is even. 

In a similar way, one can prove that eigenvalues corresponding to Hq = (2m+l)7r are 
present in the spectrum of a graph with ^-couplings on the loop, while for 5-couplings 
this is true provided the loop consists of an even number of the edges. 

If there are several halflines attached to the loop and all the b/s are equal to — 1 
or +1, respectively, we obtain the same results as before. One can easily check that for 
U = aJ + bl all entries of the energy-dependent part 

(l-A;)[/ 2 [(l-A;)C/ 4 -(A; + l)/]- 1 [/3 



of the effective coupling matrix U(k) given by (11) are equal, hence the matrix U can 
be written using multiples of the matrices J and I and the coefficient b is not energy 
dependent, i.e. U = a(k)J + bl. Since the coefficients aj can be eliminated from the 
final condition, we obtain the same results as in the energy-independent case. 
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Figure 3. A loop with two leads 

Notice that the case bj = — 1 also includes an array of edges with rationally related 
lengths and Dirichlet condition at the both array endpoints. In this case one of the 
matrices describing the coupling is Uj = diag(— 1, — 1). Similarly, bj = 1 includes the 
case of an edge array with Neumann conditions at both the endpoints, the corresponding 
matrix being Uj = diag (1,1). 

4. Examples 

As stated in the introduction our main goal is to analyze resonances which arise from 
the above discussed embedded eigenvalues if the rational relation between the graph 
edge lengths is perturbed. Let us look now at this effect in two simple examples. 

4-1. A loop with two leads 

Consider first the graph sketched in Fig. [3] consisting of two internal edges of lengths l\, I2 
and one halfline connected at each endpoint. The Hamiltonian acts as — d 2 /dx 2 on each 
link. The corresponding Hilbert space is L 2 (M + ) ©L 2 (M + ) ©L 2 ([0, Zi])©L 2 ([0, l 2 \)\ states 
of the system are described by columns ip = (gi, g 2 , fi, f2) T ■ For a greater generality, let 
us consider the following coupling conditions |ES89] which include the 5-coupling but 
allow also the attachment of the semiinfinite links to the loop to be tuned, and possibly 
to be turned off: 

A(0) =/ 2 (0), h(h) = f 2 (l 2 ) , 

A(o) =« r 1 (/((o) + ^(o)) + W 1 (o), 

h(h) = - a^(fi(h) + m)) + 7 2 ^ 2 (0) , 

<?i(0) =7i(/((0) + /^(0)) + arVi(0), 

92(0) = -7M(h) + &(h))+a; l g' 2 (0)- 

Following the construction described in Sec. [2] and parametrizing the internal edges by 
li = 1(1 — A), I2 = 1(1 + A), A G [0, 1] - which effectively means shifting one of the 
connections points around the loop as A is changing — one arrives at the final condition 
for resonances in the form 

sinfc/(l - A) sin £;/(! + A) - Ak 2 fc 1 (k) fa 1 (k) sin 2 kl 
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+ k[fc\k) + /3 2 - 1 (£;)]sin2fcZ = 0, (15) 

where 0-\k) := a" 1 + 1 « T . 

We are interested how the solutions to the above condition change with respect 
to change of the length parameter A — > X' = A + e. It is easy to check that any 
solution k depends on e continuously, and therefore for small e we can thus construct 
a perturbation expansion. Let k be solution of (15) for A and k solution for A'; the 
difference k = k — ko can be obtained using the Taylor expansion 

Kl[sm(2k l) - A sin(2A; /A)] - 4«^/? 1 ~ 1 (A;o)/£ ! " 1 (A;o) sm2k Q l - 

-^[2k^\k a )^\k,) + k 2 (fc\k )p 2 (k ) + h{kv)d 2 \k ))} sin 2 k l + 

+K{fc 1 (k ) + f3 2 \k ) + ^{k )k + k(k )k ) sm2k l + 2Klk (fc 1 (k ) + 

+(3 2 ~ 1 (ko)) cos2k l — kI[e cos k le sin k l(2X + e) + 

+ (2A + e) cos k l(2X + e) sin k le] + 0(k 2 ) = sin k l(2X + e) sin k le , (16) 

where $j(k ) = i|7j| 2 /(l — ikoaj 1 ) 2 . This equation can be used to determine k in the 
leading order. Denoting the coefficient of k by f(k ) and the rhs of the above equation 
by g(X, e) we find that the error in such an evaluation is 

j_°(« 2 )_ 1 q(9 2 {\s)\ 



f(ko) f(k ) \f 2 (k ,e) 



2\ 



Since the rhs of (16) is 0(e) as e — > 0, the error we make by neglecting the term 0(k 
is 0{e 2 ). In fact, in the vicinity of the embedded eigenvalues, i.e. for 2Xk$l close to 
= 2mr the error is even smaller, namely 0(e 4 ) as we will see below. 

In fact, we can get more from eq. (16) than just the perturbative expansion. We 
are interested in the global behaviour, i.e. trajectories of the resonance poles in the 
lower complex halfplane as A changes. To obtained them one should solve eq. (15), 
numerically since an analytic solution is available in exceptional cases only. One can, 
however, solve also numerically the approximate equation (16) starting from A = — 
where corresponding the embedded eigenvalues given by kl = nir are present, and 
taking e perturbations of the successive solutions. This method is simple and we have 
employed it in the examples below, with a sufficiently small step, e — 5 • 10 -5 . To 
check the consistency, we have compared the results in the second example with a direct 
numerical solution of eq. (15) found with the step 0.05 in the parameter A, and found 
that they give closely similar results, the relative error being of order of 10~ 3 . 

Examples of poles trajectories obtained in the described way from eq. (16) are 



shown in Figs. [4]-[6] Eq. (15) has the real solution kl = nn, nGN for A = m/n, m E N, 
the corresponding eigenfunction is ip = (0, 0, sinnirx/l, — smnnx/l) T . On Fig. [4] 
corresponding to n = 2 the pole returns to the real axis when A = 1/2 and A = 1. 
On the other hand, Fig. [5] with n = 3 shows the situation when the pole returns to the 
real axis only for A = 2/3, while for A = 1/3 and A = 1 the appropriate solution is a 
resonance. Similarly, the pole on Fig. [6] where n = 2 returns to the real axis only if 
A = 1. To show how fast the poles are moving, the change of the parameter A from 
to 1 is marked by changing the colour from red (A = 0) to blue (A = 1; visible online). 
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-0.0005 - 



-0.001 - 



-0.0015 - 



-0.002 - 



-0.0025 




6.28 6.29 6.3 6.31 6.32 6.33 6.34 

Figure 4. The trajectory of the resonance pole in the lower complex halfplane starting 
from fco = 2n corresponding to A = for I = 1 and the coefficients values a^ 1 = 1, 
a^ 1 = —2, |7i| 2 = 1, a^ 1 = 0, a^ 1 = 1, I72I 2 = 1, n = 2. The colour coding (visible 
online) shows the dependence on A changing from red (A = 0) to blue (A = 1). 







-1 - 




_2 I I I r i i 

9 9.2 9.4 9.6 9.8 10 10.2 

Figure 5. The trajectory of the resonance pole starting at ko = 3ir corresponding to 
A = for the coefficients values a^ 1 — 1, a^ 1 — 1, a^ 1 — 1, a^ 1 — 1, |7i| 2 = I72I 2 = 1, 
n = 3. The colour coding (visible online) is the same as in the previous picture. 
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Figure 6. The trajectory of the resonance pole starting at fco = 2tt corresponding 



to A = for the coefficients values a 1 1 = 1, a 2 1 = 1, o^ 1 = 1. 
I T2 1 2 = 1, n = 2. The colour coding is the same as above. 



1, l7i| 2 = 1, 



Let us now investigate the asymptotic behaviour of the resonances in the vicinity 
of the embedded eigenvalue, in particular, the angle ip between the pole trajectory 
emerging from ko = nir/l with Ao = m/n,m6 {0,1, ... ,n} and the real axis. For small 
k the difference e = A — Ao is also small. We use a rewritten form of the condition (15), 

f(k, A) = cos 2kl\ - cos2kl - 8k 2 fa 1 (k) fa 1 (k) sin 2 kl 

+ 2k(/3^(k) + fc\k))sin2kl = 0. (17) 

The function f(k, A) is, with the exception of points k - 
partial derivative with respect to A is at Aq is equal to zero, hence 



iotj, continuous and its first 



o = /(M)~M),A ) + 



d 2 f 



df 

dk 
d 2 f 



Ann 



d\ 2 



(feo,Ao) 



(fco.Ao) 



d\ 2 



fc ,A 



8k 



k ,\ 



- 4(H) 2 cos 2kl\ = -4(7rn) 



For small k we obtain using ( 16 ) 



Tin 



tany? 



Im k 
Re k 



fcol7l| 2 , 



fcp |T2 I 2 



1 + fc Q»2 



(I, 



a 9 — 



fc 2 |7i! 2 



fc 2 |72| 2 a 2 ~ 1 ' 
1+kla.- 2 



kn 



nix 
T 



(18) 



For |7i| = |7 2 | =0 the poles are real and ip = 0; this is the case when the loop and 
the leads are decopupled and the eigenvalues remain embedded. On the other hand, if 
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a. 



a 1 1 



a 



-i 



a 



then the real part of k is zero and the pole trajectory goes 
from k perpendicular to the horizontal line, i.e. (p = n/2. 

Furthermore, let us investigate the behavior of the pole trajectories hight in the 
spectrum, i.e. for large values of n. Suppose that k = k + k, k = mr/l, \k\ <C ir/l; 
then 

cos 2kl\ — cos 2kl = cos 2kolX cos 2kIX — sin 2kol\ sin 2kIX — cos 2kI 
= (cos 2nnX - 1) - sin (2?mA) 2kIX + 0(k 2 ) . 



The condition (17) for small n becomes 



TOT r 

(cos 2nnX — 1) — sin (2vmA) 2nlX + 2— [f3^\k ) + f3 2 \k ) 



2kI + OU 2 



0. 



Using the expressions of coefficients (3j(k) we obtain 

17/ 



for 5 • 1 7^ . 



^| 7i | 2 + C(l) 



for oL- 1 







The quantities appearing above, 

| cos (2n7rA) — 1| < 2 and | sin (27tkZA)| < 1 
are bounded, thus for af 1 7^ and 1 7^ 0, we have 



while for a, 1 



llm/d < 



and a 9 1 



/ 


7i 


2 /«r 2 


+ l72| 2 /«2 2 




2(nn) 2 (ar 1 + c^ 1 


- bilV^r 1 




7 2 | 2 /«2 1 ) 2 



+ C»(n" 



llm/d < 



I 



the inequality reads 
' 0(n~ 3 ) 



and for a 1 



2(nn) 2 |-yi 1 2 + I72I 2 
0, ct2 1 7^ we have 
/ 1 



Im/d < 



0{n 



2(7m) 2 | 7l | 2 

Let us summarize the discussion of this example. The poles of the resolvent are 
given by the condition ( 15 ), or equivalently, by ( 17). If A = m/n, m G N, real eigenvalues 
corresponding to kl = rnr, nsN, occur. They may correspond to a particular pole of 
the resolvent returning to the real axis for A = m/n, m E N, as in Fig. |4j However, for 
other coupling conditions, the pole may return only for certain A — see Figs. [5] and [6j 
while for other rational A its place may be taken by the pole which has been a resonance 
for A = 0. The angle between the resonance trajectory and the real axis does not depend 
on A and is given by (18). If the pole trajectory is near the original eigenvalue, then the 
distance from the real axis is of order of 0(n~ 2 ) for large n. 
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Figure 7. A cross-shaped resonator 



4-2. A cross- shaped graph 

Let us now consider another simple graph, this time consisting of two leads and two 
internal edges attached to the leads at one point - cf. Fig. [7j the lengths of the internal 
edges are l\ = 1(1 — A) and I2 = 1(1 + A). The Hamiltonian acts again as — d 2 /dx 2 on the 
corresponding Hilbert space L 2 (IR + ) © L 2 (R + ) © L 2 ([0,Zi]) © L 2 ([0,Z 2 ]), and the states 
are described by columns ip = (gi,g2, fi, f2) T - This time we restrict ourselves to the 5 
coupling as the boundary condition at the vertex and we consider Dirichlet conditions 
at the loose ends, i.e. 

/i(0) = / 2 (0) = <&(()) = <y 2 (0), 
/i(Ji) =/ 2 (/ 2 ) = o, 

«/i(0) = /((0) + ^(0)+^(0)+^(0). 

Using the same technique as above we arrive at two equivalent forms of the condition 
for resonances, k sin 2kl + (a — 2ik) sin kl(l — A) sin kl(l + A) = or 

2k sin 2kl + (a - 2ik) (cos 2kl\ - cos 2kl) = . (19) 

Let us ask when the solution is real. Leaving out the trivial case k = we get from the 
last equation two conditions referring to the vanishing of the real and imaginary parts 
of the Ihs, 

T17T 

sin2/c/ = =^ kl = —, neZ, 

TiTT T17T 

= cos 2kl A — cos 2kl = cos mcX — cos rnr = 2 sin — ( 1 — A) sin — ( 1 + A) 

=^> nX = (n — 2m), m G Z . 

Hence A = 1 — 2m/n, m G No, m < n/2. If the difference k = k — k$ is small we 
obtain from (19) 

K w — 2(a — 2ik ) sin k le sin k l(2X + e)|2i[cos 2k l(X + e) — cos2/c /] 

+ (a — 2iko)2l[(X + e) sin 2k$l(X + e) — sin 2k Q l] - 2 sin 2k l - Ak l cos 2k l } . (20) 

Similarly as in the previous example, the error here is 0(k 2 ), i.e. 0(e 2 ), and for 2Xk^l 
close to = 2n7i it is even smaller, namely 0(e A ). In the latter case, the above expression 
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6.2 6.4 6.6 6.8 7 7.2 7.4 7.6 

Figure 8. The trajectory of the resonance pole starting at ko — 2tt for the coefficients 
values a = 10, n = 2. The colour coding (visible online) is the same as in the previous 
figures. 




Figure 9. The trajectory of the resonance pole for the coefficients values a = 1, n = 2. 
The colour coding is the same as above. 
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Figure 10. The trajectories of two resonance poles for the coefficients values 
a = 2.596, n = 2. We can see an avoided resonance crossing - the former eigenvalue 
"travelling from the left to the right" interchanges with the former resonance "travelling 
the other way" and ending up as an embedded eigenvalue. The colour coding is the 
same as above. 



for ko = mc/l, A = m/n and small e yields 

2(a - 2ik ){k lef 



Ak l 



nne 



a — 2i 



nii 



I 



The slope of the pole trajectory at its start from fc is equal to 



tan <p 



Im k 2mr 



Re k 



al 



cp = arctan 



2nir 
al 



(21) 



As we have said, the embedded eigenvalues occur in accordance with (19) at 
kl = nir/2, nfZ for A = 1— 2m/n, m e No, m < n/2. The geometric perturbation gives 



rise to pole trajectories which can be found from (19), or from (20) with a sufficiently 



small step. Examples worked out using the second method are on Figs. [8-10 We see 
that a resolvent pole may return to the same point, or it may become another eigenvalue 
or a resonance. Another interesting type of behaviour, an avoided resonance crossing, 



can be seen on Fig. 10 



5. The general case 



After analyzing the above two examples, let us look what could be said about the 
geometric perturbation problem in the general case. 
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5.1. Multiplicity of the eigenvalues 

Suppose that k is an eigenvalue of multiplicity d embedded in a continuous spectrum 
of H. First we will assume that knlo — 27rm. Our aim is now to determine whether 
ko is still eigenvalue (and what is its multiplicity) if the lengths of the graph edges are 
perturbed. We will write the lengths as l'j = lo(rij + Ej) assuming that rij G N for 
j G {1, . . . , n}, while nj is not an integer for j £ {n + 1,..., iV}. 

From the construction described in the proof of theorem |3.1| we find that the 



condition (13) is not affected by small lengths variations of the "nointeger" edges, 



j G {n + 1, ... , N}. Hence the number of rationality- related eigenvalues of the perturbed 

graph referring to the first n edges does not depend on perturbations of the other 

edge lengths. The spectral condition ( 12 ) can be written as det J(k) = if we put 

J(k) := C{k) + S(k). Using the expansion 

., kUnj+Ej) . kUnj + Ej) 

ik cos — =F sin — 

2 2 

k Q l nj j k l Ej . k l Ej\ , s 

= cos — ^ I ik cos =F sin I + 0{k - k ) , 

and an analogous one for cos Ho< "j +e j) _|_ fa gj n fc; o(%+£j) Qne g nc } s the entries of J(k) 
can be rewritten as 



, . , . kolonj ( koloSj knlnEj 
Ji,2j-i{k) = {Ui, 2 j-i ~ Ui >2 j) cos — - — ik cos — sin 



+ 



ts r \ kohnj ( .. k l Q Ej . k l £j\ . . 

V>i,2j-i ~ Oifij) cos — - — Itko cos — h sin — - — I + 0(k - k ) 

T / r \ / \ k l nj ( knlnEj . k l £j \ 

Ji,2j{k) = {u i>2 j-i + Uiflj) cos — - — cos — h ik sin — - — + 



{8i,2j-\ + 5 ij2 j) cos 



k l nj 



— cos 



knlnEj 



+ ikn sin 



knlnEj 



+ O(k-k ) 



For small enough e/s and a real nonzero noninteger ko the terms cos 



k l rij 



ikn COS 



k l Sj 



sin 



k l £j 



and cos fc °'° £j + zfc sin /g °' o£j are nonzero. After dividing the columns of J(A;) by 



k l e 



these terms and using the arguments from the proof of Theorem |3.1| one arrives at the 
following conclusion. 

Theorem 5.1. In the setting of Theorem 3.1\ suppose that the rank of M even is smaller 
than 2n. Let us vary the edge lengths, l'j = ln{nj + Ej) with sufficiently small Ej 's; then 
the multiplicity of the eigenvalues e = k^ = 4m 2 7r 2 //Q due to rationality of the first n 
edges is given by the difference between 2n and the rank of the matrix 



/ un + 
«2i - 1 + e\ 
«31 
U41 



U 2JV-1,1 



«12 - 1 + eJ 

U22 + ef 

«32 
W42 



U 2JV-1,2 
U2N.2 



«13 
"23 
"33 + S| 
«43 — 1 + e£ 



U 2]V-1,3 
«2JV,3 



"14 
"24 



«44 + £9 



U 2]V-1,4 



Ml,2n-1 
«2,2n-l 
«3,2n-l 
«4,2n-l 



u 2N-l,2n-l 
U2N,2n-l 



Ul,2n \ 
M2,2n 
«3,2n 
"4,2n 



W2JV-l,2n 



«2JV,2n 



where 



(1 — fcg) sinfcoZo£j 
2ifco cos koloEj — (1 + fc§) sin fco'o^j 



gb/^\ — 2ifc o(-l + cos fc Zoej) - (1 + ^o) sinfc o'o £ J 

3 2ifco COS kglgEj — (1 + fcp) sin kgloEj 
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In a similar way one can treat the case when fco^o is equal to odd multiples of it. 
Then we employ the expansion 



ik cos 



klo(rij + Ej) 



=F sin 



kl (nj + Ej) 



sm 



k l nj 



-iko sin ^°^ £jl =F cos ^°^ £jl ) -|- C>(A; 



klojnj+ej) | ^ g - n fcZo(nj+ej) 



2 1 ™ 2 ' ' - / '-" ) 
with the help of them we arrive 



and an analogous expression for cos 
at the following conclusion. 

Theorem 5.2. In the setting of Theorem 3.2 suppose that the rank o/M dd is smaller than 
2n. Passing to Z'- = lo(nj + Ej) with small enough Ej's, the multiplicity of the eigenvalues 

e = /eg = (2m + l) 2 7r 2 /Zo due to rationality of the first n edges is given by the difference 
between 2n and rank of a matrix 



( 



Al 



«ii + f 

«21 + 1 - 
U31 
"41 



"2JV-1,1 
"2JV.1 



«12 + 1 - e 
«22 + ef 

"32 
«42 



«2JV-1,2 
U 2N,2 



"13 
"23 
"33 + 

U43 + 1 — e£ 



M2AT-1.3 
U 2JV,3 



"14 
"24 
+ 1- 
M44 + £2 



"2JV-1.4 
U 2JV,4 



«2,2n-l 
"3,2n-l 
"4,2n-l 



«2AT-l,2n-l 
u 2N,2n-l 



Ml,2n \ 
"2,2n 
«3,2n 
«4, 2n 



«2jV-l,2n , 
u 2N,2n I 



with Ej and e^ defined in previous theorem. 



5.2. Total number of poles of the resolvent after perturbation 

In general an embedded eigenvalue can split under the geometric perturbations considered here, 
a part of it being preserved with a lower multiplicity while the rest is turned into resonance(s). 
Above we have shown what the reduced multiplicity of the embedded eigenvalue is, now 
we complement this result by showing that the total number of poles produced in this way, 
multiplicity taken into account, remains locally preserved. Before stating the result, let us 
first demonstrate two useful lemmata. 

Lemma 5.1. Let (k,E) 1— > g(k,e) : C x IR m — > C be a function uniformly continuous in e 
for all e G U £o (0) and k G Wr(&o), £0 > 0, R > 0, and holomorphic in k in Un(ko) for all 
£ G U £o (0). Furthermore, let lim^o g(k,s) = {k — ko) d . Then there exist such 5 > and 
£q > that for all £ G U e i (0) the sum of the multiplicities of zeros of g(k, s) in Us(ko) is d. 

Proof. Since g is holomorphic, we have the Taylor expansion 

00 

g (k, e) = J2 a p( £ )( k ~ k o) P = P(k, £) + {k- k ) d+1 h(k, £) = P(k, £) [1 + (k - ko)h(k, £)] , 

p=0 

where P(k, e) is a polynom of order d in the variable k, furthermore, lim^o h(k, e) = 
and limg-^o h(k, e) = lim^o^ — ko) d h(k,£)/P(k,£) = 0. Due to the fundamental theorem 
of algebra P(k, e) has d zeros, not necessarily different, whose distance from ko depends 
continuously on e. On the other hand, we have \/53e' : Ve G U £ i Q (0),\/k G Wr(/so) : \h(k,£~)\ < 5 
in view of the above limit relations; choosing then 5 < 1 / R we can conclude that zeros of the 
term [1 + (k — ko)h(k, £)] lie outside the ball 1An(ko). □ 

The following lemma slightly generalizes the result to a larger class of g(k,s). 

Lemma 5.2. Let (k,£) i— > F(k,£) : C x IR m — > C be a function uniformly continuous in e 
for all £ G U £o (0) and k G Un(ko), £q > 0, R > 0, and holomorphic in k in Uji(ko) for all 
£ G U £o (0). Suppose that F(k,0) has in Wij(fco) o single zero of multiplicity d at the point ko; 
then there exist such 5 > and e' > that for all e G U £ i^ (0) the sum of the multiplicities of 
zeros of F(k,£) inUs(ko) is equal to d. 
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Proof. In view of the holomorphy of F and the fact that F has a zero of order d in ko one 
has F(k, e) = (k — ko) d f{k,e), where lim^o f(k, e) ^ 0. Because / is continuous in e we have 
f(k, e) ^ for all e* G W £ ^(0), k € Upi{ko). Hence / does not contribute to zeros of F in IAn(ko) 
and Lemma 15. II can be used. □ 



This conclusion allows us to demonstrate the indicated result. Our aim is to determine 
the number of resolvent poles, multiplicity counting, of the quantum graph with perturbed 
edge lengths in the neighbourhood of an original pole of multiplicity d. In particular, we want 
to find out whether the number of solutions of the condition ^ — into which we substitute 



from (11) — changes in the neighbourhood of fep- I n the notation of the previous lemma, the 



function F is given by the Ihs of ^ and the vector e describes the change of the edge lengths. 

Theorem 5.3. Let Y be a quantum graph with N finite edges of the lengths li, M infinite 

edges, and the coupling described by the matrix U = ^ j^ 1 j^ 2 J , where U4 corresponds to the 

coupling between the infinite edges. Let ko satisfy det [(1 — ko)Ut± — (1 + ko)I] 7^ and let ko 
be a pole of the resolvent (H — Aid) -1 of a multiplicity d. Let T £ be a geometrically perturbed 
quantum graph with the edges of lengths /j(l + e) and the same coupling as V . Then there 
exists an Eo > such that for all e £ U £Q (0) the sum of multiplicities of the resolvent poles in 
a sufficiently small neighbourhood of ko is d. 



Proof. One can rewrite the condition ( 12 ) for poles of the resolvent into the form F(k, e) = 0, 
where e is the vector of differences of the lengths of the internal edges. Using the form of the 
matrices D\(k) and D2(k) and Eq. ( fIT| ) one can easily check that if det [(1 — ko)!!^ — (1+fcoK] 7^ 
then there exists a neighbourhood Un{ko) where F(ko,e) is holomorphic in k and uniformly 
continuous in e, hence Lemma |5.2| can be applied. □ 



Notice that the condition det [(1 — ^0)^/4 — (1 + &oK] 7^ is automatically satisfied for 
ko G IR + because of the inequality \(ko + l)/(&o — 1 ) | > 1 and the fact that the eigenvalues of 
C/4 do not exceed one in modulus. 



Acknowledgments 

The research was supported by the Czech Ministry of Education, Youth and Sports within the 
project LC06002. We thank the referee for suggestions which helped to improve the text. 



References 

[AC71] J. Aguilar, J.-M. Combes, A class of analytic perturbations for one-body Schrodinger operators, 

Commun. Math. Phys. 22 (1971), 269-279. 
[EL06] P. Exner, J. Lipovsky, Equivalence of resolvent and scattering resonances on quantum graphs, 

Adventures in Mathematical Physics (Proceedings, Cergy-Pontoise 2006), 447, Providence, R.I., 

(2007), 73-81. 

[ES89] P. Exner, P. Seba, Free quantum motion on a branching graph, Rep. Math. Phys. 28 (1989), 
7-26. 

[ET07] P. Exner, O. Turek, Approximations of singular vertex couplings in quantum graphs, Rev. Math. 

Phys. 19 (2007), 571-606. 
[Ku04,05] P. Kuchment, Quantum graphs: I. Some basic structures, II. Some spectral properties of 

quantum and combinatorial graphs, Waves in Random Media 14 (2004), S107-S128; J. Phys. A: 

Math. Gen. 38 (2005), 269-279. 



Perturbed quantum graphs with rationally related edges 



24 



[Ku08] P. Kuchmcnt, Quantum graphs: an introduction and a brief survey, Analysis on Graphs and its 
Applications, Proc. Symp. Pure. Math., AMS 2008, pp.291-314. 



